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ABSTRACT 

Blazar spectral models generally have numerous unconstrained parameters, leading to ambiguous 
values for physical properties like Doppler factor 5-d or fluid magnetic field B' . To help remedy this 
problem, a few modifications of the standard leptonic blazar jet scenario are considered. First, a log- 
parabola function for the electron distribution is used. Second, analytic expressions relating energy 
loss and kinematics to blazar luminosity and variability, written in terms of equipartition parameters, 
imply (5d, B' , and the principal electron Lorentz factor 7' The external radiation field in a blazar 
is approximated by Lya radiation from the broad line region (BLR) and sa 0.1 eV infrared radiation 
from a dusty torus. When used to model 3C 279 SEDs from 2008 and 2009 reported by Hayashida et 
al. (2012), we derive i5d ~ 20 - 30, B' ~ few G, and total (IR + BLR) external radiation field energy 
densities u ^ 10"^ - 10""^ erg cm""^, implying an origin of the 7-ray emission site in 3C 279 at the 
outer edges of the BLR. This is consistent with the 7-ray emission site being located at a distance 



R < V'^ctyar ^ 0.1(r/30)^(i^ar/10 s) pc from the black hole powering 3C 279's jets, where t^jar is 
the variability time scale of the radiation in the source frame. For low hadronic content, absolute jet 
powers of « 10% of the Eddington luminosity are calculated. 
Subject headings: Gamma rays: general ~ Galaxies: jets - radiation mechanisms: nonthermal 



1. INTRODUCTION 

The time-varying emission spectrum of a blazar en- 
codes important information about particle acceleration 
and radiation, jet structure and environment, and the 
mechanisms by which black holes power jets. Determin- 
ing the jet physics and environment from the spectral en- 
ergy distributions (SEDs) of blazars is a tricky problem 
of inversion. Part of the problem is that observational 
limitations make it difficult to obtain a detailed picture 
of the blazar SED, as the blazar continuum extends over 
some 20 orders of magnitude, from MHz radio frequen- 
cies to TeV 7-ray energies, and therefore requires coor- 
dinated campaigns at many wavelengths. The observing 
program furthermore has to be as nearly simultaneous as 
possible, due to the highly variable nature of blazars. In 
spite of these obstacles, extraordinary progress has been 
made to obtain a relatively complete description of the 
blazar SED. 

The EGRET instrument on the Compton Gamma Ray 
Observatory first showed that 100 MeV - GeV emission 
often dominates the apparent power output of blazars 
during fiaring states of flat spectrum radio quasars (FS- 
RQs), and that the 7 rays form a spectral compo- 
nent that is distinct from the nonthermal synchrotron 
emission radiated at radio, IR and optical frequencies 
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(see. e.g.. lH"artman et al']|1992l: iWehrle et al.lll998[ ). The 
ground-based Cherenkov 7-ray telescopes demonstrated 
that the BL Lac class are powerful sources of TeV 
emission (IPunch et al.lll992[). and tha t some FSRQs — 
3C 279 dM AGIC Collaboration' '2008), PKS 1510-089 
( Wagner et al 2010: Cortina 2012 ), and 4C -t-21.35 



Aleksic et al .1120 111) — are even detected at very-high en- 



ergies (VHE; > 100 GeV). 

The Large Area Telescope (LAT) on Fermi, which 
nominally operates in scanning mode, has revolution- 
ized our understanding of blazars by providing unin- 
terrupted light curves and detailed SEDs of blazars. 
In campaigns involving radio, IR, optical. X-ray, and 
ground-based 7-ray telescopes, superbly detailed SEDs 
have been measured that challenge our theoretical under- 
standing. A few examples include the high-syn chrotron- 
peaked BL Lacs Mrk 421 (jAbdo et al.ll2011bD and Mrk 
501 (Abd o et al. 2011cf) . the intermediat e synch rotron- 
peaked BL Lac 3C 66A (lAbdo et all l2011dD . and 
the FSRQs 3C 454 .3 (|Abdo et al.l I2009D and 3C 279 
(jAbdo et all l2010aD . Mrk 421 and Mrk 501 are 
reasonably well fit with a leptonic synchrotron/SSC 
model, whereas an additional Compton-scattering com- 
ponent is required in FSRQs and low- and intermediate 
synchrotron-peaked BL Lac objects. 

A standard blazar paradigm consisting of a mag- 
netized plasma ejected with relativistic speeds in a 
collimated outfiow along the polar axes of a rotating 
black hole has been developed to exp lain the overall 
obser vational features of blazars (see IBoettcher et all 
I2OI2I for a recent detailed review). But even the simple 
synchrotron/SSC model suffers from a proliferation 
of parameters, which greatly increases with the intro- 
duction of external radiation fields having uncertain 
spectrum and energy density. The range of possible 
fits to the SEDs makes it difficult to extract unique 



and therefore meaningful blazar jet properties. Part of 
this difficulty has to do with characterizing the non- 
thermal electron spectrum, and another part with the 
range of possibilities of external p hoton sources, which 
can include the accretion-disk (iDermer et al. 19921 

isen 



7' ~ 7p therefore contain the maximum nonthermal elec- 
tron energy per unit logarithmic interval A (log 7') in the 
electron distribution. Thus 



Dermer &: Schlickeiseil 119931 : IDermer fc S chli ckeii 

2002|) , broad lin e region (BLR; [S ikora et al. 199 4), and 



mole cular torus (iBlazeiowski et al. 2000; Arbeiter et al.l 
2OO2I: [S okolov & Marscher 2005) radiations (see also 



Bottchei1l2007l : iSikora et al.ll2009l) . 



To help alleviate these problems, we adopt a specific 
model where the target radiation fields are the BLR ra- 
diation and IR radiation from a dusty torus surrounding 
the nucleus. We employ a log-parabola function in Sec- 
tion 2 to describe the electron energy distribution in the 
comoving fluid frame. Equipartition relations between 
the energy densities of the magnetic field, the nonther- 
mal electrons, and the observed properties of the blazar 
synchrotron spectrum are derived in Section 3. In Sec- 
tion 4, model results are shown for target photons of 
the external radiation described by monochromatic fre- 
quencies and various energy densities. Spectral fitting of 
recent SEDs of 3C 279 (Hayashida ct al. 2012) are used 
to illustrate the method. 

In a companion paper ([Cerruti et al.l [2013D . we im- 
prove the model for the external radiation field, replac- 
ing monochromatic infrared emission by thermal spectra, 
and the BLR radiation field by a complex of atomic emis- 
sion lines rather than by the Lya line alone. We show 
that scattering of the BLR radiation under equipartition 
conditions can explain the GeV spectral cutoff discov- 
ered with Fermi-LAT in 3C 454.3 and other FSRQs and 
low- and intermediate-synchro tron-peaked BL Lac ob- 
jects (|Abdo et al.ll2009ll2010bl ). 

2. LOG-PARABOLA ELECTRON DISTRIBUTION 

In leptonic blazar jet models, the synchrotron and 
Compton 7-ray emission components are radiated by 
nonthermal electron distributions that are assumed to be 
isotropic in the jet fiuid frame. One technique is to fit the 
data by injecting power-law particle distributions and al- 
lowing the particles to evolve ii i response to radiative and 
adiabatic energy loss es (e..s... JBottcher fc Chiand 120021 : 
iModerski et all 120031 ). This approach requires many pa- 
rameters for the cutoff energies, injection indices, and 
power, but is potentially useful to follow dynamic blazar 
spectral behavior. 

Contrary to this approach, we abandon any preconcep- 
tions about particle acceleration, and employ the sim- 
plest functional form that is able to provide reasonably 
good fits to the snapshot SED data. [3 For this purpose, 
we assume that the log-parabola function 

logW^N^ij')] = log[7;2iVe(7;)] - b log(777;)' (1) 

gives a good description of the nonthermal electron 
Lorentz- factor (7') distribution N^{'-f') in the comoving 
fiuid frame, where the electron distribution is assumed to 
be isotropic. The entire description of the electron en- 
ergy distribution is then given by three parameters: the 
normalization, the peak (or principal) Lorentz factor 7' 
and the spectral curvature parameter b. Electrons with 
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^ IFinke et al.l 1)20081 ) uses a similar approach for SSC modeling 
of TeV blazars. 



that is, the electron distribution function is a smoothly 
curving function with continuously varying spectral in- 
dex. Normalizing to the total comoving nonthermal 
electron energy £'^ through £'^ — ra^c^ j^ d-f' j' N^{'-f') 

gives S'^ ^ mec^j'p^ Ki%)Iib), where /(&) = ^/^TlnlO/b 

(VTrlnlO ~ 2.69). Note that the relationship is not ex- 
act, because we have replaced 1/7!, with in the lower 
limit of the integral. 

The log-p arabola function ha s been used to de- 
scribe blazar (iMassaro et al.ll2004[ ) and gamma-ray burst 
(jMassaro et al.ll2010D spectra. It has also been recog- 
nized as a convenie nt functional form f or the electron 
energy distribution ("Massaro et aL| [2006( ). and is treated 
in detail byiTramaccrc et al. (20111). 

3. EQUIPARTITION RELATIONS 

One goal is to determine, by comparison with blazar 
SED data, whether the synchrotron and 7-ray spectra 
are reasonably well fit by a log-parabola function for 
the electron energy distribution, within the scenario that 
the blazar jet is described by a standard relativistic jet 
model with quasi-isotropic electrons trapped in a ran- 
domly oriented magnetic field with mean strength B' in 
the fluid frame. The principal comoving electron Lorentz 
factor 7' functions like the minimum Lorentz factor in 
a fast-cooling blast- wave scenario (jSari et al.lll998| ). but 
here only represents instantaneous or time-averaged (not 
evolving) particle and photon distributions. 

3.1. Parameters and Derived Quantities 

In the stationary frame of the black-hole/jet system, 
we can list the important observables that describe the 
blazar SED: 

1. Ls ~ 10*^L48 erg s~^, the apparent isotropic bolo- 
metric synchrotron luminosity; 

2. Eg = hvs/meC^ = 8.13 x 10~^i/i4, the peak syn- 
chrotron frequency, that is, the frequency at which 
the vlv synchrotron spectrum reaches its maxi- 
mum value, in units of meC^, where 10^*z/i4 Hz 
is the peak synchrotron frequency in the source 
frame; 

3. tyar — 10*^^4 s, the variability time scale, where 
tvar = i°v^ 1^ + -z), and i°^* is the measured vari- 
ability time scale; 

4. Aec = Lec I Ls^ the external Compton-dominance 
factor giving the ratio of the apparent isotropic 7- 
ray luminosity, assumed to be from Compton scat- 
tering of external isotropic radiation, to Lg ; and 

5. Assc = Lssc/Ls, the internal, or self-Compton- 
dominance factor, giving the ratio of the SSC lumi- 
nosity to Ls- This ratio is more difficult to deter- 
mine from the blazar SED, because the SSC com- 
ponent is concealed by the EC component, though 
it may be reveal at hard X-ray energies. 



A complete SED contains additional information in the 
detailed features of the spectrum. 

Within the framework we set up, the derived syn- 
chrotron/SSC/EC model parameters are: 

1. B' , the comoving magnetic field, in G; 

2. (5d, the Doppler factor; 

3. 7' the principal electron Lorentz factor; 

4. r^ — cdj^tyan the comoving radius of the emission 
region; 

5. Mo, the external radiation field energy density; and 

6. eo = hvo/nieC^ , the dimensionless external radia- 
tion photon frequency. 

The dynamical crossing time associated with the 
Schwarzschild radius of a IO^Mq black hole is sa 10^ 
s. For only a few extremely bright blazar flares does 
the Fermi L AT have sensitivity to probe to a few 
hour time scale (lAbdo et al.ll2011at iTanaka et all 120111: 
iFoschini et al.l 1201 1[ ), and detecti on of much more rapid 
variability from TeV blazars (e.g.. lAharonian et al.ll2007l : 
lAleksic et al.l[201ll ) represents a major puzzle. To keep 
the exposition uncluttered (but incomplete), we let t^ = 
1 in our calculations, keeping in mind that this parameter 
should be determined by observations of the variability 
properties of the blazar at the time that the observing 
campaign took place. 

The comoving magnetic field energy density u'g, = 
B''^/8tt. In the standard blazar model, we can normal- 
ize the total energy density uj^j of particles, field, and 
radiation in the outflowing plasma blob through the ex- 
pression 

""tot — "S'll + Cpar + Cph) = 



^B'y^ + Ce + Cp/nuc + Cs + C*) 



(3) 



The subscripts here refer to magnetic field (-B')j includ- 
ing both large-scale and turbulence and MHD waves; 
particles (par), divided into nuclear (p/nuc) and elec- 
tronic (e) components; and photons (ph), consisting of 
the internal synchrotron (s) and external (*) radiation- 
field energy densities in the comoving frame. In the 
Thomson regime, Assc = Lssc/Ls = u'Ju'b' ~ Cs, 
and Aec = Lec/Ls = u'^/u'^, « C*- 

The simplest equipartition relation is Wg = Ce^s'? Ce — 
1, which also minimizes jet power. We let Cp/nuc — 1 
in the following calculations, which is the upper limit 
to a jet with low baryon content. The baryonic content 
affects only the total power and not spectral properties; 
see Section 6.4. 

3.2. Constraints and Equipartition Relations 

One advantage of using the log-parabola function, be- 
sides having only three parameters to characterize the 
nonthermal electron distribution, is that in the limit 
h -^ 00, the electrons become monoenergetic. Here we 
solve in that limit, and use the analytic results to provide 
parameters as input for a numerical model. 



A first constraint arises from simple kinematics. In a 
blob geometry, the comoving synchrotron energy density 
u'g = L^™/47rrj^c(5Q, where r^ = cS-otvar, so 

Lf- = 47rrf c<5^M^ . (4) 

A second constraint, which is exact in the limit 6 — >■ 00, 
is the synchrotron luminosity constraint 

Lr = ^caTMi,,7Vo7;f<S^. (5) 

When the equipartition relation rrieC^NQ^' /Vl = u'^ = 
(eu'g,, with Vj^ = Anr'^/S, is substituted into the syn- 
chrotron constraint, eq. ([5]), and the result is equated 
with the luminosity given by the kinematic constraint, 
eq. dU, we find 



^'^r^'^'^'^^^^'r' 



p 



(6) 



In the 6 — >■ cx) asymptote, the peak synchrotron fre- 
quency can be expressed in terms of the critical magnetic 
field Bcr = m1c^/eh as 



'^ = 2 BTr^^ '^° 

from which, with eq. (jH]), we derive 



Sr,B'^ 
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Replacing this result in eq. (0]) gives 
2i3 



(7) 
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Inverting eq. ([7]) to get an expression for B' as a function 
of fo, and using eq. ([9]) for (5d gives 
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Using eqs. ©, ([S]), and pHt . we have 

5/8.1/4 



^48 
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(10) 



(11) 



3.3. EG Component 

In the Thomson regime, the peak frequency of the ex- 
ternal Compton-scattered component is 

^■EC _ 4 , ivi ^ -- o 9 x^ in3 Ce 3/2^1/4 n- 

^C - 3 \P^lp) ^* = ^-^ X lU -^Ja ^li ^48 V M ELya ■ 

(12) 
The onset of the Klein-Nishina correction for external 
Compton scattering takes place when 47' e' > 1. Tak- 
ing e' = (5De*, the photon energy in the stationary frame 
where KN effects become important is at ec,KN = l/12e*. 
For Lya photons with e* = 2 x 10~^, -Ec.kn — 2.2 
GeV, and this could be the rea son for the spectral cut- 
off in FSRQs and LSP blazars (jAckermann et al.|[2010l: 
ICerruti et al1[20Tl . 



3.4. SSC Component 

The peak SSC Compton 7-ray frequency is approxi- 
mately given by 
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The onset of Klein- Nishina effects in SSC scattering takes 
place at photon energy es > es.KN = S^^/'^lp, implying 
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(14) 



(15) 



Thus, KN effects are most important in high-synchrotron 
peaked (HSP; vf"" > W^ Hz) BL Lac objects. 

3.5. Synchrotron Self- Ah sorption 

We use the (5-function approximation to calculate syn- 
chrotron self-absorption (SSA). The SSA coefficient for a 
photon of frequency m^c^e' /h in the fluid frame is given 
in this approximation by 



-TT Acre 
36 e' 



7s 






Is 



e'/2e', 



B ' 



^ , (16) 

(jDermer fc MenonI [20091 eq. (7.144), with corrections). 
Here Ac = h/rUeC is the Compton wavelength, r^ = 
e'^/nieC^ is the classical electron radius, and e'^ — 
B'/Bcr- Taking ^(7') = N'^{i)/V{,, and substituting 
Equation ^ for N'Xl) gives 



SL 



Ac re 



18 mec2/(6)Fj,'7;4 



(2 + 61ogx)x" 



-(4+6 log x) 



(17) 
where x = ^e'/2e^/7p. 

The SSA opacity r^/ — 2Kg/r^. The unabsorbcd 
spectrum is multiplied by the factor m(t) = 1/2 + 
exp(— T)/r — [1 — exp(— t)]/t^, with t — t^i, to give the 
absorbed spectrum. 

3.6. Jet Power 
We calculate absolute jet power Ljet using the relation 



(18) 

for a two-sided jet (jCelotti &: FabianI 119931: 
iCelotti fc Ghisellinil I2008D . where the total internal 
energy density uj^j is given by Equation ([3]). This 
equation is surely valid for isotropic particles and 
randomly oriented fields, but the photon distributions in 
the comoving frame have peculiar isotropics not taken 
into account by Equation ()18p , which make the different 
angular distrib utions of syn chrotron and Compton 
radiation fluxes (|Dermeiill995f) . 

For photon powers, we therefore use expressions re- 
lating apparent isotropic and absolute photon fluxes 



(jDermer et al.|[20T2l . Appendix A). The apparent {Liso,s) 
and absolute {Labs,s) synchrotron fluxes are related by 
the expression 



^abs.s — 81 Liso,s/ iO-Q , 



(19) 



and a similar expression likewise holds for the SSC fluxes. 
For external Compton (C) fluxes, the relationship is 



^abs,C 



32r iiso,c/5(5D I 



and 



Lph — Labs.s + Labs.C ■ 



(20) 



(21) 



These very different behaviors are important for power 
requirements in photon-dominated jets. The kinetic 
power 2T£'^/tyar required by the jet to transport the ra- 
diating electrons and magnetic field is found to be much 
smaller than the photon power, so that baryon-loading 
docs not have a big impact on jet power requirements 
until it is > 10 x the electron energy. 

4. NUMERICAL RESULTS 

Using the derived values of (5d, B' , 7^, and r^, along 
with the log-parabola description of the electron spec- 
trum. Equation ((ij, we numerically calculate vL^ SEDs. 
Th e methods f or cal culati ng the SEDs a re pre sented 
in iFinke et all (poM ) and iDermer et aP (pOOl . We 
use the head-on approximation for Compton scatter- 
ing (jjonesl I1968D . which is valid for relativistic elec- 
trons (7 ^ 1) scattering lower-energy photons (e ^ 
7), with Klein-Nis hina effects fully taken into account 
(jDermer fc MenonI 120091 ). 

In the model calculations shown in this section, we cal- 
culate SEDs for external isotropic photon sources with 
e» = 2 X 10~^ characterizing the BLR (10.2 eV), and 
e, = 2 X 10^^ (0.1 eV) describing radiation from the 
dust torus. The former value is motivated by the strong 
Lyg line in blazars, specifically, 3C 454.3 ([Bonnoli et al.l 
120111) . and the latter value by calculation s of AGN 
emissi ons reprocessed by a dust clouds Nen kova et al.l 
(2008,?). The peak power of the reprocessed emission oc- 
curs at w 10/i (the 10/i silicate feature), with broad wings 
from « 3/i to ~ 30/z. This peak frequency correspond to 
e, « 2 X 10^^, or an effective blackbody temperature of 
w 440 K, characteristic of warm dust. 

Relevant for characterizing the IR environment of 
blazars are r ecent observations of th e unusual FSRQ 4C 
-F21.35 by. iMalmrose et al.l (|2011| ). They analyzed 5 
- Sbfj. Spitzer, Sloan Digital Sky Survey, Two Micron 
AU Sky Survey, and Swift UVOT data of 4C -^21.35, 
decomposing its spectrum into a nonthermal power law 
and two-temperature dust model, made of a warm dust 
component with effective temperature T ~ 660K radi- 
ating w lO''^ erg s~i, and a hot dust component with 
T « 1200 K radiating « 8 x 10"*^ erg s-\ both on the 
pc scale. We examine hot dust fits for 3C 279 below, 
but consider warm dust (e* = 2 x 10^^) for an initial 
parameter study. 

When fitting actual blazar spectra, both IR dust and 
BLR photon sources contribute and have to be appropri- 
ately adjusted to give a good fit. Detailed fitting requires 
a range of atomic lin es and multiple thermal components 
(jCerruti et al.|[20T3l ). but this two- line simplification il- 
lustrates this approach to modeling the SEDs of blazars 
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" ^e = Cep = 1; blank entries in Tables are equal to first vertical 

entry above. 

' Equipartition synchrotron/SSC model 



TABLE 2 
Blazar model output 
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log [Ljet(iph) 








(G) 




(erg 
cm-^) 


(10^5 
cm) 


(erg s-i)] 


la 


17.8 


17.8 


6.75 


141.2 


0.0431 


5.33 


46.6(47.1) 


lb 


8.89 








0.172 




46.0(46.5) 


Ic 














46.0(46.8) 


Id 


17.8 








0.0431 




46.6(47.4) 


2a 


11.5 


11.5 


7.80 


163. 


0.136 


3.46 


46.0(46.5) 


2b 


5.77 








0.545 




45.4(45.9) 


2c 














45.4(46.1) 


2d 


11.5 








0.136 




46.6(47.1) 


3a 


7.49 


7.49 


9.00 


1.8802 


0.431 


2.25 


45.4(46.1) 


3b 


9.99 


9.99 


3.80 


7.94e2 


0.0431 


3.00 


45.1(45.9) 


3c 


13.3 


13.3 


1.60 


3.35e3 


0.00431 


4.00 


44.9(45.4) 


4a 


6.49 


6.49 


4.38 


9.17e2 


_ 


1.95 


43.9(44.4) 


4b 


8.65 


8.65 


1.85 


3.97e3 




2.59 


43.7(44.0) 


4c 


11.5 


11.5 


0.78 


1.63e4 




3.46 


45.4(43.6) 


4d 


15.4 


15.4 


0.33 


6.88e4 




4.61 


43.2(43.3) 


4e 


20.5 


20.5 


0.14 


2.90e5 




6.15 


42.9(43.0) 



in order to derive basic properties of the jet and environ- 
ment. 

Figure [T] shows the resuhing model SEDs for values 
characteristic of a very powerful FSRQ blazar with a 
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Fig. 1. — Model blazar SEDs from the radio to the 7-ray im- 
plied by near-equipartition conditions, using synchrotron luminosi- 
ties and peak frequencies characteristic of a powerful FSRQ like 3C 
454.3. In this calculation, fs = 0.5 and f» = 10. The bolometric 
apparent synchrotron luminosity Ls = 10** erg 3- '^. The jets are 
viewed off axis [6 = 1/F in a and d, and on-axis in b and c. The 
target photons are Lya (e* = 2 X 10—5) in a and b, and IR dust 
photons (e* = 2 X 10—'') in c and d. Separate synchrotron, SSC, 
and EC components are shown for a. Thin solid line shows the 
synchrotron emission component in the absence of self-absorption. 



bolometric synchrotron luminosity Ls = 10^* erg s~^ 
and peak photon frequency Vg = 10^'^ Hz. Input values 
for this and subsequent figures in this section are given 
in Table 1, and derived values are listed in Table 2. In 
the calculations, we use the relation 



1 



(1 + N^)Su 



(22) 



to relate F and (5d, where N-p — T6, valid in the limit 
T ^ 1, 9 <^ 1. We consider two cases: iVp = 1, cases 
[T^ and [TJi, and A^r = 0, cases [TJd and [TJ;. We also use 
a single external monochromatic line appropriate to the 
BLR (e» = 10.2/511000) in cases a and b and an IR torus 
(e* = 2 X 10~^) in cases c and d. 

Already we see an interesting result: a break at GeV 
energies is formed when jet electrons scatter Lya/BLR 
photons, but not when scattering IR photons, which in- 
stead make vF^ peaks near 10 MeV. If the electrons in 
the flaring jet are within the BLR, they scatter both 
IR and BLR photons, whereas when the radiating jet 
is found far outside the BLR, scattering of IR photons 
dominates. So an FSRQ with a GeV break would also be 
accompanied by a 7-ray external Compton dust feature. 
In contrast, flares occurring far outside the BLR would 
make a blazar SED peaked at MeV energies, as seen in 
some low- frequency peaked blazars, e.g., CTA 102 and 
PKS 0528+134, as indi cated by joint OSSE, COMP TEL, 
and EGRET analysis (jMcNaron-Brown et al.lll995D . 

The different beaming factors of the synchrotron and 
SSC components on the one hand and the EC compo- 
nents on the other are seen when the i/L^ SED for 9 = 
1/T,Sb =ris compared with the SEP for 9 = 0,5b^ 2F 
(|Dermedll995HGeorganopoulos et al.ll200lD . In Figure [U 
F « 18 for 6* = 1/F and F « 9 for 9 = 0, with photon 
energy densities uq w 0.043 erg cm-'^ and uq = 0.17 erg 
cm"-^ corresponding to SSC and EC factors Cs = 0.5 and 
C* = 10, respectively. A larger energy density is required 




16 20 

logv{Hz) 



Fig. 2. — Same as for Fig. [T] except that Lsy„ = 10* 
characteristic of powerful blazars like 3C 279. 



for the on-axis jet model, because the imphed F factor is 
smaUer. 

Even though (s = 0.5 in Figure [U the SSC luminos- 
ity is about equal to the synchrotron luminosity. This 
is because the log-parabola function is not a simple 5- 
function, and the electrons with 7' > 7' are more effec- 
tive scatterers than lower energy electrons, insofar as the 
Thomson energy- loss rate |7t| oc 7'^. This also accounts 
for shifting of the model peak frequency to i^s > 10^^ Hz. 
In the limit that b — >■ 00, i/g = 10^^ Hz. For the same rea- 
son, the peak frequency of the Compton-scattered dust 
component gets shifted to higher energies than expected 
from simple Thomson estimates. When 9 w 1/F, the 
Thomson-scattered 7-ray luminosity L^ of the dust com- 
ponent is about C*(= 10) X larger than Ls, but is even 
larger when viewed down the jet axis due to the different 
beaming factors. The Compton-scattered Lya compo- 
nent has a smaller apparent luminosity than that of the 
dust component because of the onset of Klein-Nishina ef- 
fects. The absolute jet luminosities for these cases almost 
reach the Eddington luminosity for a 10^ Mq black hole 
when 9 = 1/F, but are significantly less for the on-axis 
viewing. 

Figure [2] shows a calculation similar to Figure 1 , but 
with an order-of-magnitude smaller apparent bolomet- 
ric synchrotron luminosity, now with Lg — 10^^ erg s^^. 
The values of F become smaller than in the cases shown 
in Figure 1 , but the energy densities of the external radi- 
ation fields must become larger in order to produce the 
same Compton dominance (see Table[2] The absolute jet 
powers become smaller, though not as fast as decreasing 
(X ij^. Again, the equipartition solution (Ce = 1) with 
ti — 1 gives a breaking GeV spectrum from target Lya 
photons, while the upscattered IR photons peak at ~ 1 - 
10 MeV in a lyF^ SED, while breaking at slightly higher 
energies. IR radiation from hot dust would increase the 
break energy by a factor ~ 3, still well below the GeV 
range. 

Figure [3] shows calculations for parameters character- 
istic of I ntermediate Synchr otron Peaked (ISP) blazars, 
defined (|Abdo et al.l l2010c| ) as blazars with observed 
peak synchrotron frequencies in the range 10^^ Hz < 
1/°'"* < 10^'^ Hz. The synchrotron and external Comp- 
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Fig. 3. — Near-equipartition models with input parameters 
characteristic of ISP blazars like 3C 66A and W Comae. Here 
Lsyn = 10*^ erg s~^, fs = 0.5 and (J, = 10, and a range of val- 
ues of Us = 0.1, 1, and 10 are chosen. Vertical lines show the range 
bounding peak synchrotron frequencies Us defining ISP blazars, de- 
fined by 10^* Rz < Us < 10^^ Hz, defined according to lAbdo et al.l 
pOlOd ). Thick lines are for Lya photons, and thin lines are for 
1200 K IR photons from the dusty torus. Separate components, 
including the unabsorbed synchrotron emission, are shown for case 
a with target Lya photons. 



ton parameters are kept as before, with (g = 0.5 and 
C* = 10, but the value of vu is allowed to vary from 
0.1 to 10. The different SEDs for target IR and Lya 
radiation fields are shown by the thin and thick curves, 
respectively. Notice that values of Vg that give 10^^ Hz 
< i^s < 10^^ Hz are really about an order of magnitude 
smaller for 6 = 1, as explained above. The discrepancy 
decreases with increasing 5, and becomes exact in the 
limit 6 — > oo. 

Depending on the source redshift, the cases shown in 
Figure [3] could all qualify as ISP blazars. The increase 
in Vg coincides with an increase in F, but with a reduc- 
tion of Mo because of the increased efficacy of Compton 
scattering an external isotropic radiation field when F 
is larger. The onset of Klein-Nishina effects reduces the 
importance of the Compton-scattered Lya field, so that 
the SED is dominated by 7 rays produced from target 
IR torus photons. In one of the cases shown, the peak 
and cutoff from scattered IR photons is near 1 GeV. If 
in fact, the dust torus was the dominant target photon 
source for the blazar 7-ray SED, then a linear correla- 
tion between the 7-ray and synchrotron peak frequencies 
is predicted. 

For these ISP-type parameters, structure may be found 
in the GeV regime due to the competition between SSC 
and Compton-scattered Lya photons. There is also the 
possibility of an inverted spectrum when the EIC emis- 
sion starts to dominate the SSC emission. Such a fea- 
ture in ay be present in the SED of 3C 66A (jAbdo et al.1 
l2011dD . but will require better data and fitting to estab- 
lish. 

Figure 2] shows equipartition models for blazars with 
characteristics similar to those of HSP BL Lac objects. 
Comparison of the SEDs of M rk 421 and Mrk 501 from 
multiwavele i igth ca mpaigns of lAbdo et al.l ()2011bD and 
'Abdo et al.l (|2011c[ ). respectively, with models |4}: and|4}l 
suggest that these TeV BL Lac objects will be well fit 
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Cases 4 a-e, 
roughly adapted 
to Bl Lac objects 




TABLE 3 
Blazar input parameters for 3C 279" 



Epoch 


Lis 


un 


e/ij 


(s 


CLya 


Or 


Ld^sk iW^" 

ergs 1) 


A 


0.06 


0.01 


2 X 10-^ 


0.06 


8 


4 


0.55 


A" 






5.4 X 10-^ 




6 


12 




B 


0.15 


0.02 


2 X 10-'^ 


0.04 


14 


11 


1.0 


B" 


0.1 


0.05 


5.4 X 10-7 


0.1 


18 


8 


0.5 


C^ 


0.18 


0.02 


2 X 10-'^ 


0.04 


7 


2.2 


0.2 


J)C,d 


0.18 


0.03 




0.045 


16 


5 


0.5 


"ti = Ce 


= (p/nuc = Nr = b = 1; EBLR = 2 


X 10-5 






''Hot dust solutions: T = 1200 K, em = 5.4 
'^No hot dust solutions with t4 = Q = 1 
''In this fit, b = 1.1 


X 10-7 







16 20 

logv(Hz) 

Fig. 4. — Near-equipartition SSC models characteristic of nearby 
BL Lac objects like Mrk 421 and Mrk 501. The model shows results 



for Ls 



10* 



Cb = 0.5, and 



;i^i4 



1,2,3,4,5, corre- 



sponding to curves a through e, respectively. Vertical lines define 
the range of average peak synchrotron and Compton frequencies 
for Mrk 421 and Mrk 501, and light and dark arrows point to 
values of u^, vLi,^, vq, and uL,^^ for Mrk 4 21 and the Mrk 501 , 
resp ectively, from the c ampaigns reported in lAbdo et al.l Il2011bl ) 
and lAbdo et all 11201 Id l . 



with our equipartition model. The broader widths of the 
synchrotron and Compton SED components in Mrk 501 
compared to Mrk 421 resulting from a smaller b in the 
log-parabola description of the electron energy distribu- 
tion can help explain the differences in vq and uL^q 
between these two sources. 

For this parameter set, we also calculated the energy 
density of Lya radiation needed to make a significant 
contribution to the 7-ray fluxes. When ULya ^ 10~* erg 
cm ' ' ' ... 



normalized such that /^ deL{e) = Ldisk, so that 
r(4/3) = 1/1.12. The i/F^ spectrum of the accretion-disk 
is therefore /°„1, = eLdjsfe(e)/47rd|, where e = (1 + z)e°'"' 
and dL is the luminosity distance. The value of tmax 
depends on the spin of the black hole and relative Ed- 
dington luminosity, but for simplicity we let emax = 10 
eV, typical of the characteristic temperature of the UV 
bump in Seyfert galaxies. 

The spectral luminosity of the IR dust component is 
approximated by a thermal distribution normalized to 
the IR luminosity Ljf(. Thus 



£LiR{e) 



IbLiR (e/ey 



7r4 exp(e/e) - 1 



(24) 



with corresponding vFi, spectrum f^J^^ = eLj]^(e)/4:TTdj^. 
Here 6 = k^T/meC^ = 2 x IQ-^, taking T = 1200 
K. A dust covering factor of 20% is assumed, so that 
Lib. = 0.2Ldisk- The energy density of the IR radiation 
is limited by the energy density of a blackbody, namely 



the emission would make an EC component with 



l^C 



> 10 



41.5 



erg 



just barely visible in the i^Li, 



Ubbierg cm"3)(r) = 9.36 x lO^^e^ 



3x 10-^( 



T 



440 K' 



0.016( 



T 



1200 K' 



(25) 



spectrum. The Lya energy densities in Mrk 421 and 

Mrk 5 01 are, however, far below this level. iStocke et al.l 

(|2011| ) report Lya luminosities of 2.4 and 5.2 xlO^° erg 

s-i for Mrk 421 and Mrk 501, respectively. The total Thus the energy density at most reaches w 10"^ erg 

jet power is « lO^'* erg s"^ for the models in Fig. H and cm ^ for warm dust, and w 0.02 erg cm ^ for hot dust. 

the accretion-disk luminosity Ld should be greater than 

this if most of the jet energy comes from mass accretion 

(black-hole rotational energy could supply a comparable 

amount). Even for a BLR as small as 0.01 pc, the Lya 

energy density ULyn ~ 10^°L4n/4:TTRn m„^c « 3 x 10^^ 



5. MODELING THE SEDS OF 3C 279 



'^ya 



erg cm~ , far smaller than what is necessary to make 
the EC 7-ray flux brighter than the SSC flux. 

The low Lya luminosities of Mrk 421, Mrk 50 1, and 
also PKS 2005-489 and PKS 2155-304 ( StockeeTall 
I2011D . two other HSP BL Lac objects, validates the pre- 
sumption that TeV blazars have a very tenuous BLR 
environment. By contrast, FSRQs can h ave very large 
Lya luminosities; in the case of 3C 454.3 (|Bonnoli et al.l 
1201 ID . reaching si (2 - 4)xl0'^5 g^g s'^. 

A complete model SED requires, for consistency, the 
emission spectrum from the accretion-disk and from the 
dust. The accretion-disk spectral luminosity is assumed 
to be described by a Shakura-Sunyaev spectrum given by 



f-Ldisk{e) = 1.12X, 



disk 



4/3 



exp(-e/emaa;) , (23) 



iHavashida et al.l ()2012D have organized 3C 279 cam- 
paigns around Fermi-LAT with wonderful results. The 
data in Figures [5] and [6] show SEDs from quasi- 
simultaneous observing campaigns for four periods of 
Fermi-LAT observation, namely Epochs A (MJD 54682 

- 54729; 4 Aug 2008 - 19 Sept 2008), B (MJD 54789 

- 54809; 19 Nov 2008 - 9 Dec 2008), C (MJD 54827- 
54877; 27 Dec 2008 ~ 15 Feb 2009), and D (MJD 54880- 
54885; 18 Feb 2009 - 23 Feb 2009). The SEDs com- 
prise X-ray data from Suzaku|!| Swift XRT, XMM New- 
ton and RXTE, optical/UV data from Kanata, GROND 
and Swift UVOT (170-650 nm), IR data from Spitzer, 
and radio data from CARMA and OVRO. These SEDs 
are contemporaneous and not strictly simultaneous, as 
different telescopes require different observing intervals 
to achieve the displayed sensitivity. Nevertheless, these 

^ Suzaku consists of the XIS (0.3-12 keV), HXD/PIN (10-700 
keV), and HXD/GSO (40-600 kcV). 
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Fig. 5. — Near-equipartition blazar model fits (by eye) to the 
SEDs of 3C 279 (Havashida ct al. 2012) for Epochs A (upper panel) 
and B (lower panel), with input parameters given in Table [3] and 
implied quantities from the model in Table HI Dot-dashed curves 
and components, peaking from low to high-frequencies, are the 
optically thin and absorbed nonthermal synchrotron emission, the 
thermal dust emission, the accretion-disk radiation, and SSC, EC 
dust, and EC Lya emission components assuming a warm dust 
target with ejn = 2 X 10~^. Dashed curves show modeling with 
hot dust, tjn = 5.4 X 10"'^. Insets show detail of fits at X-ray 
energies. 



SEDS are some of the most detailed FSRQ spectra yet 
measured. 

Figures [5] and H] also show by-eye fits to 3C 279 from 
our equipartition model. This approach is restricted to 
sources with known redshift, and the re dshift of 3C 279 
is pa rticularly well known (z = 0.5362; iMarziani et al.l 
I1996D . Putting in observables t^ar, Lgyn and Vpk,syn 
gives 5y>, B' and r^ through the equipartition formula- 
tion, Equations (|9]), ([T0|. and (|TT|) . respectively. With 
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Fi g. 6. — N ear-equipartition blazar model fits to the SEDs of 3C 
279 IIHavashi da et al. 2012) for Epochs C (upper) and D (lower), 
with input parameters given in Table [3l and derived quantities in 
Table \^ Components as in Fig. \E\ 



the goal of limiting the number of parameters, we take 
i4 = Ce = Cp/nuc ^ Nr = b ^ 1, and vary only L48, 

^14, Cs, Cec and Cdust (Cec + Qust = (*)■ In fact, tyar 
can be measured, as can Nr if the angle to the jet axis 
can be determined, and such values should be used when 
available. The value Cp/nuc only affects the jet power, 
not the spectrum, as is proved by examining Equations 
Q, dUl), and dni - dni)- Relaxing Ce away from unity, 
which increases the jet power, can be used as a last re- 
sort, but for the purposes of the study of a system in 
electron/magnetic field equipartition, we keep Ce = 1- 

In the fitting, there is some freedom to assign the 
accretion-disk temperature and luminosity, and there- 
fore dust luminosity for the assumed 20% covering factor. 



TABLE 4 
Model output values for 3C 279 



Epoch 


r,<5D 


B' 


7; 


"Lya 


UIR 


log[Ljet(iph) 






(G) 




(10-3 erg 
cm-3) 


(10-3 erg 
cm-3)/ 

T„,n(K)" 


(erg s-l)] 



A 

A* 
B 
B* 
C= 



20.7 3.2 60.4 



30.8 
21.4 
31.8 
31.8 



1.8 
2.7 
1.8 
1.7 



93.4 
143 
92.3 
116 



5.6 

4.2 
1.4 
8.8 
0.64 
1.3 



2.8/776 
8.5/1024 
1.1/610 
3.9/840 
0.22/410 
0.41/480 



46.3(46.0) 
46.4(46.2) 
46.7(46.3) 
46.5(46.4) 
46.4(45.9) 
46.7(46.2) 



"Minimum blackbody temperature to exceed u/^j 

''Hot dust solutions 

'^No hot dust solutions with t4 = Ce = 1 

Sometinies the accretion-disk emission is important to fit 
structure in the optical, but in other cases gives only an 
upper limit to disk luminosity Ldisk- 

The results of fitting, by eye, the data for Epochs A-D 
of 3C 279 are shown in Figures [S] and [5] and listed in Ta- 
bles [3] and m The spectral components and total spectra 
are shown in these figures for the case of a warm dust IR 
field {eiR = 2 x 10^^). The different components, as writ- 
ten in the figure caption, are the unabsorbed and syn- 
chrotron self-absorbed nonthermal synchrotron spectra, 
the IR graybody field, the Shakura-Sunyaev accretion- 
disk field (big blue bump), the SSC, EC 7 rays from dust 
and, at the highest energies, EC from Lya. For compar- 
ison, we also show fits in Fig. [S] with a hot-dust external 
IR target radiation field (e/_R — 5.4 x 10^). In both these 
cases, we were able to find hot dust solutions, facilitated 
in the Epoch B case by its poorly defined synchrotron 
component. 

We could not find solutions that gave acceptable fits 
to Fig. inifor the hot dust solution, provided we keep the 
other parameters, in particular, ^4 and ^e fixed to unity. 
In these cases, the Compton-scattered IR component had 
little effect on the rising X-ray component, and the SSC 
radiation peaked at too low of a frequency to match the 
Swift XRT data. The model fits with the warm dust 
appear to simulate well enough the blazar spectra. In 
Epochs A and B, we can also obtain good fits with hot 
dust. In fact, the dust emission spectrum is far broader 
than a black body, so more detailed examination must 
consider scattering a broad-band IR spectrum. 

6. DISCUSSION 

In the Fermi era, the blazar SEP can be detailed 
enough as in the case of 3C 279 (jHavashida et al.ll2012D 
that each model component is required and constrained 
by the data: the nonthermal synchrotron component, of 
course, for the nonthermal radio/mm/IR/optical radia- 
tions; accretion-disk radiations to explain the optical/UV 
blue bump; nonthermal SSC X-rays; hard X-rays ob- 
served by the Swift XRT and BAT and NuStar made by 
SSC and Compton-scattered IR (EC-IR) photons origi- 
nating from outside the jet; and GeV 7 rays by external 
Compton scattering of Lya and BLR photons (EC-BLR) . 
If the models lacked the EC-IR 7-ray component, the 
model SEDs would be very different, both at the lower 
end of the Fermi-LAT range and in the hard X-ray/soft 
7-ray spectrum explored by Swift BAT, Suzaku, INTE- 



GRAL, and NuStar. Note that accretion-disk radiations 
are required to obtain good fits to the optical spectra in 
Epochs A and B, and the model does not apply to the 
< 10^^ Hz radio measurements. 

6.1. Model Implications 

In this approach to multiwavelength FSRQ modeling, 
fitting to a blazar SED implies values of (5d, B' , 7^, and 
r^ for a specific value of tyar and observing angle 9. Our 
fits to the multi-epoch data from 3C 279, Figs. [5] and 
m imply r « 5d ~ 20 - 30 and B' « few G, assum- 
_ ing 9 — 1/T and ^4 = 1. These values are not too dif- 
ferent fr om earlier determinatioii s of th ese numbers in 
work by iBlazeiowski et ahl ()2000r n and iHartman et al.l 
(|2001| n when fitting to multiwavelength data associated 
with earlier EGRET campaigns, where reasonable fits 
were obtained by guessing values for F and B' . 

Assuming that the technique implying F and B' give 
results that correspond to the actual conditions of the 
radiating jetted plasma expelled by the black hole, then 
the energy densities of the surrounding radiation fields 
are derived with an uncertainty refiecting the range of 
parameters that yield a good fit to the data. It is not, 
however, easy to say what that uncertainty is without a 
detailed numerical study. Here our aim is more limited: 
to demonstrate the ability of the equipartition model to 
extract jet and enviromental parameters. 

6.2. Location of ^-ray production site 

Epochs A-D range in duration from 6 days (Epoch 
D) to 3 weeks (Epoch B) and 6 weeks (Epochs A and C), 
and statistics accumulated during these observing times 
show up in the err or bars of the Fermi-LA T SEDs shown 
in Figs. [5] and H (jHavashida et al.ll2012D . The derived 
target photon energy densities for Epochs A - D of 3C 
279 are MLya ~ 10^^ - 10^^ erg cm"^ for the Lya field 
and u/fl ~ 10~^ - 10^** erg cm"'^ for the IR field, ranging 
over a factor of ~ 40 for the different fits. The total 
energy density of the external radiation field is « 10^^ 
- 10"^ erg cm"'^. The Lya rad iation comprises n early 
one-third of the broad line flux (jTelfer et al.ll2002t l. and 
the total broad line flux may amount to « 10% of the 
accretion-disk luminosity reprocessed by BLR gas. 

Simple scaling relations from reverberation mapping 
studies indicate that the BLR radius i?BLR can be 
written in terms of the accretion-disk luminosity Ldisk 
through the expression 



^BLR — O.I1 



L, 



disk 



1046 



erg s 



- pc 



(26) 



(|Ghisellini fc Tavecchioll2008D . If a fraction t < 1 of the 
radiation is scattered or reprocessed into line radiation, 
then Equation ((26|) implies that the energy density of the 
scattered radiation is 

I erg cm . (27) 



^scatt 



4^-R|LRf 



0.031 



^0.1' 



so that the typical energy density of Lya radiation in the 
BLR is at the level of (0.01 - 0.1) erg cm"'^, consistent 



"^ B' = 0.81 G, r = 20, 6* : 
8 B' = 1.5 G, r = 5 - 



l/E: 

15, t 



7br = 150 
observing angle from 



ILahteenmaki fc ValtaoTal I I1999I) 



10 



with the derived energy densities of Epoch A. At later 
Epochs, the inferred Lya energy density is well below 
w 0.01 erg cm"'^. 

If the emission region is at the outer edges of the BLR, 
in this case at > 0.1 pc, then the derived Lya energy den- 
sities could be consistent with all Epochs A-D, including 
the hot-dust solutions for Epochs A and B. Note that 
the IR energy densities for a characteristic temperature 
are constrained by Equation ((25|) . For the derived IR 
energy density, we give the minimum blackbody temper- 
ature, Tmim in Table m From the derived energy den- 
sities shown in Table HI we see that the largest required 
IR energy densities for a hot dust solution is compatible 
with the blackbody temperture, whereas the warm dust 
solutions are marginally compatible with the blackbody 
limit. The warm dust solutions in Epochs C and D are 
compatible with the blackbody limits. 

For the hot dust solution, the energy density of 1200 
K blackbody radiation is sa 0.02 erg cm~^, from Equa- 
tion ([25]) . The radius of a 1200 K blackbody radiating 



10''"'^L45 erg s^^ of thermal radiation is Rbb ~ Q.li\/L^ 
pc. The smaller derived values of the IR dust energy den- 
sity, uiR = 10~^u_3 erg cm^-^ (0.2 < u_3 < 6), are con- 
sistent with the same luminosity, but from a larger region 
of characteristic size Rir = \/ubb/uiB,Rbb ~ O-^/y/uZ^ 
pc. 

The IR luminosity that we assume for 3C 279 is small 
by comparison with other sources. For 4C +21.35, an 
FSRQ like 3C 279, thermal IR radiation with a luminos- 
ity of 7.9 ± 0.2 X 10^^ erg s~^ was measured, as was a 
weak detection of dust from the FSR Q CTA 102 with 
IR lu minosity of 7 ± 2 x 10**^ erg s~"^ (Malmros e et al.l 
120111) . On the other hand, an upper limit on dust emis- 
sion of 2.3 X 10*5 erg s^^ was found for PKS 1510- 
089, another FSRQ. The IR luminosity derived by fit- 
ting the 3C 279 SED (Table IS, 2./^ - 1 - 2xl0'*5 
erg s~^, is somewhat lower than these values, but this 
number could be increased by factors of several by as- 
suming a smaller dust covering factor, without affect- 
ing the model fits. We can write Wdus* [erg cm^^] = 
3 X 10-^{LiR/10'^^erg/s)/{R/10^^ cm)^, with the char- 
acteristic size scale of about 1 It-yr for the IR in 3C 279. 

Energy densities of target radiation fields can vary 
from one epoch to another for two reasons: (i) by chang- 
ing locations of the emitting jet plasma to places with 
different ambient radiation fields, or (ii) by having the 
jet 7-ray production take place at a similar location while 
the surrounding radiation environment changes. With a 
dynamical time of ~ 1 yr for the IR field, and variations 
in ujn by > 10 from the warm dust solutions over the 
seven months comprising the 4 Epochs studied, case (i) 
seems preferred. Indeed, we can use Epoch A, for which 
the largest value of ujr, m_3 « 6 was derived, to argue 
that the characteristic size of the IR torus could be as 
small as ~ 0.2 pc. The 7-ray emission at later Epochs 
could be made at > 0.5 pc from the black hole, where 
the IR field is smaller, consistent with the smaller energy 
densities of radiations inferred in later Epochs. 

6.3. Jet Model 

Within the colliding shell paradigm, the collision ra- 
dius Rcoii ^ 2rlctvar ^ 0.2(r/30)^t4 pc, which is also 
at the outer perimeter of the BLR (Eq. (|26|) ). Collision 



radii of relativistic plasma ejecta seem consistent with a 
location at the edge of BLR, which is also on the ~ 0.1 
- 0.5 pc size scale of the IR dust emission where um 
changes significantly. In Epoch A, when the derived F 
is smallest, so Rcoii < ^Ictyar ^ 0.03(r/17)^t4 pc, the 
probable location of the emission region is deep within 
the BLR and IR torus. If the IR torus is on size scales as 
small as 0.2 pc, which depends o n dust properti es, then 
the size scales are similar. From iNenkova et al.l ([2008), 
the inner radius due to dust sublimation is 



1500 K' 



(28) 



where the bolometric AGN luminosity is 10 erg 
and T^iih i s the dust sublimation temperature. 



iMalmrose et al.l ()2011| ) argue for dust sublimation ra- 
dius of 1 - 2 pc. It could also be that the IR zone is 
not much larger than the BLR, as might follow natu- 
rally if the torus and BLR clouds are the same accretion 
flow found at different distances from the nuclear contin- 
uum, with properties defined by whether the accreting 
flow was inside or o utside the dust-sublimation radius 
(jNenkovaet al.ll2008D . 

6.4. Jet Power 

The absolute 2-sided jet power is given by Equation 
IHl), L,,t = 2^rf/3r2cuJ„t, with uj^^ = ^'^,(1 + Ce + 
Cp/nuc + Cs + C*) fo'' the total energy density, from Equa- 
tion ([3]). The jet power ijet(erg s~^) is reported in Table 
H for the 4 Epochs of 3C 279. We also give the photon 
power Lph = Lph^syn + Lph^EC, from Equations (fT9l) and 
(f20|) for the synchrotron/SSC and EC photon powers, re- 
spectively. A value of Cv/nuc = 1 is assumed in the values 
of Ljet given in Table |41 The maximum jet power found 
in our calculations, again with baryon-loading unity, is 



^jet 



1046.7 gj.g g-i « 5 X 10"^ erg s'\ 



This number should be compared with th e Edding- 
ton luminosity. From emi ssion line studies, iGu et aP 
(|2001h : IWoo fc Urr^ (|2002f l infer a black hole of mass 
(3 - 8) X IO^Mq in the 3C 279 system. The implied Ed- 
dington luminosities are in the range (4 - 10)xlO*^ erg 
s~^. It is worth remarking that if Eddington conditions 
limit maximum jet power, then the most powerful flare 
is at the Eddington limit, provided Cp/nuc is sufficiently 
small. The absolute jet powers for the warm-dust solu- 
tions, in units of 10^^ erg s~^, the maximum Eddington 
luminosity, are Cp/nuc = (1, 10, 100) = (0.14, 0.2, 0.95) for 
Epoch A, (0.54,0.71,2.5) for Epoch B, (0.26,0.46,2.4) 
for Epoch C, and (0.46, 0.63, 2.3) for Epoch D. To keep 
the absolute jet luminosity < 10^^ erg s^^ in all cases 
requires a baryon-loading factor Cp/nuc ^ 30. The pho- 
ton power is the dominant component of the the power 
budget in the jets of 3C 279; see Table H 

The apparent isotropic bolometric 7-ray luminosity ex- 
ceeds « 10''® erg s""' in Epoch B, whereas the absolute jet 
power, dominated by photons, is ~ 5 x 10^^ erg s""', for 
low hadronic content, Cp/nuc — 1- Reduction of power 
requirements is a consequence of the jet opening angle 
6 K, l/F, which leads to beaming factors /& ^0.1% for 
20 ^ r < 30. The reduction from apparent to absolute 
jet power should follow the solid angle reduction accord- 
ing to the beaming factor /{, = 1.7 x 10"'^ in Epoch A, 
and 9 w 1/30 rad; fb^6x 10^^ in later epochs. The 
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reduction in power is not by the full beaming factor be- 
cause it depends on the beaming cone of the radiation. 

6.5. Particle Acceleration and the Highest Energy 7 rays 

As a fitting function, the log-parabola spectrum seems 
to be capable of fitting multiwavelength blazar data rea- 
sonably well, even with b fixed to unity. Such curv- 
ing functions are well known to arise in second-order, 
stochastic acceleration s cenarios (e.g.. Park & PetrosianI 
119951: iBecker et al] |20061. If the quality of fits are better 
and if fewer parameters are used in the current frame- 
work than employed in shock/injection scenarios, we may 
conclude that the curving spectrum is a better model for 
the data. 

A shock framework may still be required in a second- 
order acceleration scenario, because turbulence induced 
downstream of a shock might be required for acceler- 
ation. A turbulent scenario might also naturally arise 
in some version of models involving jet reconnection 
in compact kn ots far (> pc) from the nucleus (e.g., 
iGiannios et al.l 2009), and Poynting-dominated jet mod- 
els. An interesting question is why the 6=1 solution is 
a good first approximation to the electron spectrum. 

With regard to the fits to 3C 279, it is clear that 
in Epochs A, B, and C, the model underproduces the 
^ GeV emission measured by Fermi-LAT. This defi- 
ciency cannot be cor rected by usi ng a better BLR model 
(ICerruti et al.l 120 13[ ). and the IMAGIC Collaboration! 
(I2OO8D detection of VHE photons from 3C 279 seems to 



require either leptonic scenarios far out of equilibriunfl 
or hadronic processes (Bottchcr et al.l [20091 ). In view of 
emission extending to the VHE range, a new spectral 
component may be required, and the best candidate is 
one made by UHECRs accelerated in blazars. As we have 
seen, the jet powers can accommodate large baryon-load 
factors, (p/nuc £ 30, so 3C 279 or, indeed, other power- 
ful F SRQs, could d isplay high-energy hadronic emission 
tails (jAtovan fc De rmcr 2003). 



TZeB'r'^ 



6.6. Hadronic Acceleration and UHECRs 

From the IHillasI (J1984D condition, E„u 
ZecBTSutvar = 1-6 X 10^°Z{B'/2 G)(r/30)((5D/30)i4 
eV. Thus 3C 279 is not prevented from accelerating pro- 
tons to super-GZK energies by this basic requirement. 
FSRQs like 3C 279 are, however, unlikely to be the so- 
lution to the origin of the UHECRs, insofar as no FS- 
RQs are found within the GZK radius of ~ 100 - 200 
Mpc. In spite of their great power, the low space den- 
sity of FSRQs and powerful FR2 radio galaxies makes 
them unlikely to be the primary source class power- 
ing the UHECRs. For tha t, FRl radio galaxies and 
BL Lac ob jects arc favored (jDermer fc Razzaaud 120101 : 
IMurase fc" Takami 2009). These types of blazars have 
trouble, however, accelerating protons to energies > 10^^ 
eV, above which th e composition must become heavy 
(jMuraseet al.ll2012D . 

UHECR production in FSRQs like 3C 279 might be 
revealed by cascade radiation induced by photopion pro- 

® There is no difficulty in making strong TeV radiation with 
an equipartition model; see Fig. |4] Strong flaring VHE or TeV 
emission in an FSRQ would therefore need to be accompanied by 
a synchrotron radiati on compon ent peaking at X-ray energies, as 
in the Mrk 501 flare l lPian et aLlirggSl . 



duction from escaping UHECR protons. At redshift z = 
0.5362, angular size distance dA — 1301.8 Mpc {Hq ~ 71 
km (s-Mpc)~^, r^m = 0.27 = 1 — 17a). Depositing energy 
on a distance scale of 100 - 200 Mpc with a jetted beam 
pointed nearly along our line of sight will create through 
photopion and photo pair processes a s teady, e xtended 7- 
ray e mission feature ()Essev fc Kusenk o 2010; E ssev et all 
120101 ). Under certaing conditions, the UHECR beam can 
escape from the structured re gions surroun ding the 3C 
279 nucleus without dispersal (|Murasell2012[ ). 

A blazar like 3C 279, in the ecliptic plane, at lati- 
tude 0.2°, is not preferentially located for IceCube ob- 
servations, but otherwise would be a promising target 
for PeV neutrino searches. The external Lya radiation 
field provides a target for ultra-relativistic protons to 
make a neutrino spectrum down to PeV energies (see 
Fig. 3 in (iDermer et al.ll2012D : contrary to statements of 
ICholis fc Hooped (120121) '). 

6.7. Blazar Divide and Sequence 

One of the most impressive blazar correlations in 
Fermi-LAT blazar studi es is the 7-ray sp ectral index 
versus 7-ray luminosity (|Abdo et al.l 12010a ). A partic- 
ularly noticeable separation between the BL Lacs and 
FSRQs in this plane is referred to as the blazar divide 
l|Ghisellini et al.|[2009[ ) , and interpreted in terms of differ- 
ent accretion regimes, expressed as an Eddington ratio. 
The object 3C 279 occupies, no surprise, a position con- 
sistent with a soft-spectrum FSRQ, with F = 2.31 below 
the 3.5 GeV break and a ste eper spectrum, F = 2.95, 
above ([Havashida et al.|[2012[ ). 

When the blazar sequence (jFossati et al.l 119981 ) and 
also blazar divide are fit by simple blazar mod- 
els, one obtains a ty pical electron Lorentz factor (7) 
(jGhisellini et al.lll998[) . and the general relation is that 
the peak electron Lorentz factor Y j^ ^ I'p ^ il) is 
inversely correlated with Lj, where Lj is the average 
apparent isotropic 7-ray power) , though the detaile d be- 
havior is unclear and in dispute (jMever et al.ll201l"l ). 

Obtaining an inverse correlation of 7' and Lj is not 
obvious in a first-order acceleration scenario involving a 
wind of power L = L^ > Lj, where maximum particle 
energy ex L^'^. In cooling scenarios, the cooling Lorentz 
factors go to lower energies when the radiation field is 
more intense (" Ghisellini et al.lll998l : IBottcher fc DermeH 
2002; Finkc 201^^ 

In a second-order acceleration scenario, the inverse cor- 
relation 017' and Lj can be obtained if the scale of turbu- 
lence Xmax oc 1/fcmm grows with increasing L-y. To show 
this, note that the most important losses are synchrotron 
and Thomson, with loss rate —jsyn.T = Ksyn,TLj^, as- 
suming S^/Stt -|- L^^ cx L. The energy-gain rate from gy- 
roresonant acceleration of electrons or ions with Alfvenic 
turbulence is 



q-l 



ckLnPACirLK^.n)'-'!"'-' (29) 



(jDermer et al.|[T996l) , where q is the spectrum of the tur- 
bulence and r^ = mc^ /QB' . The dimensionless Alfven 
speed /3^ = va/c — B' j ^jAwrh^nhpC? , so that if B'^ ex L^, 
then the equilibrium Lorentz factor 7' increases with L 
unless turbulence with greater total energy is injected at 
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larger size scales (e.g., Sedov radius ex EeLp)- 

6.8. Spectral Issues 

In the models for 3C 279, synchrotron self- absorption 
invariably sets in at ~ 10^^ Hz. Even at the sub-mas 
scale, some slowly varying optically thin radiation must 
be made as the jet plasma expands and radiates. Our 
one- zone model omits the complex radio-emitting geom- 
etry composed of a superp osition of self-abso rbed zones 
and shock- in-jets fcf.. lMarscher fc Ge ar 1985). A clearer 
picture of the transition between self-absorbed and opti- 
cally thin jet plasma can be made with upcoming obser- 
vations of 3C 279 using the Sub-Millimeter Array (SMA), 
which operates at 230, 345, and 420 GHz. A regular mon- 
itoring would indicate how variability properties differ 
between these three wavebands, and would make a joint 
CARMA/OVRO/SMA analysis possible. Radio variabil- 
ity properties can be further compa red with data from 
15 GHz VLBI MOJAVE monitoring (jLister et al.ll2009[ ). 

At X-ray energies, the search for inverted spectra at X- 
ray frequencies, possibly due to the transition from dom- 
inant SSC to EC fluxes with increasing photon energy, 
should be a high priority. The spectral hardening should 

Spear not only in 3C 279, but in most bright FSRQs, 
and is required from a comparison of Swift and Fermi 
blazar spectra. The quality of the X-ray data already 
challenges the model fits, as can be seen from the insets 
in Figs. [S] and m As can also be seen from the insets, the 
inverted spectrum sets in between about 10 and 40 keV 
in the 3C 279 SEDs, which is within the Swift XRT range 
and at the lower end of the BAT energy range. Measure- 
ment of X-ray polarization will also help discriminate 
between SSC and EC processes (Krawczynski 2012). 

As to what new information NuStar brings, note that 
it is much more sensitive in the 6-80 keV range than 
Suzaku, and it is in this range where we predict a hard- 
ening from the EC-IR component taking over from the 
SSC component. Noting already the quality of the X-ray 
data, high-quality 3C 279 spectra from NuStar should 
clearly show details o f spectral emission c omponents at 
hard X-ray energies (jHarrison et al.|[2010[ ). In this re- 
gard, it is important that the limitations of the fit at 
X-ray energies seen in the details of Figs. [5] and |6] can 
be mitigated if there is an emission component from the 
hot inner disk near the black hole that shines at hard 
X-ray energies, like the Shakura-Sunyaev disk can shine 
through at optical/UV frequencies. 

At hard X-ray/7-ray energie s, the CGRO analysis 
by iMcNaron-Brown et al] (|1995f ) showed that the 7-ray 
peaks are broad, and that the vFi, peaks of these sources 
are at ^ 1 MeV energies, given the hard spectra for the 
FSRQs 3C 273, 3C 279, CTA 102, PKS 0528+134, and 
3C 45.3 at OSSE (- 50 keV - 1 MeV) energies. Within 
the context of this model, the Compton-scattered dust 
component is essential to explain flux at OSSE levels, 
which should also be detectable with INTEGRAL, and 
so would strengthen the case for photons of an IR radi- 
ation field being Compton-scattered by jet electrons to 
7-ray energies. 

Though we have explored a particular paradigm, it re- 
mains to be confirmed that the 7-ray emission from the 

10 Swift BAT spectra of 3C 454.3 and 3C 273 show apparent 
hardenings above 200 MeV l|Aiello et al.|[20T2 ). 



hard X-ray to ^ 10 MeV is scattered IR radiation. Scat- 
tered accret ion-disk radiation gives a similar emission 
comp onent (jHartman et al.|[2001l : iDermer fc Schlickeiseii 
l2002t ). but is further complicated by the precise location 
of the emitter from the accretion disk. Collisions at var- 
ious distances from the black hole change dramatically 
the relative importance of the accretion-disk component. 
One thing that is not fully understood is why there is 
such a large range, by a factor of « 10 for the warm- 
dust solutions, in energy densities um for the IR field 
implied by the fits. This can be explained if the emission 
site is on the outer edge of the IR radiation field, where 
the scattering efficiency rapidly decays. The greatest IR 
energy density is recorded in Epoch A, implying an IR 
core of radius sa 0.2 pc. This is not so different from our 
estimates for the BLR radius, tblr ^0.1 pc, and for the 
collision radius Vcoii ~ 0.2 pc. If the size of the IR region 
inferred from IR luminosity and model energy density to 
be ^ 0.2 pc, then collisions outside this radius can result 
in a large reduction of received flux, which could account 
for the wide range of IR energy densities in Table |4l 

6.9. 77 opacity constraints 

Blazar observables give, within the confines of the one- 
zone approach, a minimum bulk Lorentz factor V from 
the constraint 77 —> e . If sources are opaque to 7 ray 
escape, then pair re-injection can make a synchrotron 
or Compton feature, and can initiate further cascade 
generations. With all its apparent might in the record 
flare of November 2010, 3C 454.3 was transparent to 
its own radiation with a bulk Lorentz factor F > 14 
(jAbdo et al.ll2011al ). This shows, however, that solutions 
with r,(5D^10-15 also need to include a check on the 
level of their internal opacity. 

The large, F > 20 bulk Lorentz factors found in the 
equipartition model for 3C 279 are consistent with it 
being transparent to GeV 7 rays, so that opacity and 
cascade emission can be neglected. On the other hand, 
hadronic secondary 7 radiations can be severely absorbed 
by the IR dust field (jDermer et al.|[20T2l ). A further check 
on the model is to make sure that 7 rays can escape the 
emission region by avoiding internal 77 absorption with 
the synchrotron radiation, and external 77 absorption 
with the Lya and BLR radiation. 

7. SUMMARY 

A minimalist prescription drove us to a few-parameter 
blazar model with which to fit blazar spectra, giving 5d, 
B' , and external radiation field energy densities ujji and 
ULya for the IR torus and the BLR field, respectively, as 
outputs of the fitting procedure. Even in this approach, 
the simple blazar model has 10 parameters: (1) 6d, (2) 
F, (3) B', (4) Si, (5) 7;, (6) esLfl, (7) ublr, (8) cir, (9) 
u/fl, and (10) 6, which is simplified to 6 parameters by 
assuming 9 = 1/F, 6 = 1, €blr = f-Lya = 2 x 10~^, and 
e/fl corresponds to either 440 K dust (warm dust, ejR = 
2 X lO"'^) or 1200 K dust (hot dust, em = 5.4 x 10"^). 
The six parameters are determined by the observables 

(1) Vs, (2) I/L^^, (3) tyar, (4) VSSC, (5) I'ii/sso (6) vc, 

and (7) vL^^\ values for parameters (4) - (7) depend 
on modeling the overlapping components in the X-ray 
and 7-ray spectra. An optimum model implies energy 
densities of the target radiation fields. 



13 



From fits to 4 epochs of 3C 279, taking 6 = l/F, 
tvar = 10"* s, and Qg — 1, we derive Lya radiation-field 
energy densities ULya ^ 0.001 - 0.01 erg cm~^. This is 
below values expected within the BLR, ~ 0.01 erg cm~'^, 
implying that the 7-ray production zone is instead out- 
side the BLR radius Rblr ~ 0.1 - 0.2 pc. The wide 
range of inferred IR energy densities, differing by more 
than an order of magnitude over the different epochs in 
the warm-dust solutions, is difficult to explain if the IR 
radiation field has to vary by such a large factor over the 
< 7 month period considered here. It therefore seems 
more likely that the emission regions occur outside the 
BLR, at i? ^ 0.1 - 0.5 pc, consistent with relativistic 
shell collision radii. Gradients in the IR fields in the re- 
gion containing the BLR in 3C 279 are likewise on the 
scale of the BLR. 

This gives a picture of a system with large IR gradients 
at the sub-pc scale, comparable to the size scale of the 
BLR. This might be accomplished most effectively if the 
molecular torus and accretion disk were part of the same 
gaseous flow toward the central nucleus. 

Within the constraints of the equipartition model pre- 
sented here (assuming that ^4 = Ce = 1), warm-dust solu- 
tions were found for all four Epochs A - D, but were only 
marginally consistent with the blackbody constraint im- 
plied by Tmin in Table[3]in Epochs A and B. Hot-dust so- 
lutions were found for Epochs A and B, made possible by 
the lack of MeV constraints in Epoch A and synchrotron 
constraints in Epoch B. In Epochs C and D, only warm- 
dust solutions were possible, and hot-dust solutions were 
excluded because the target photon frequency was too 



large to make any useful spectral contribution at hard 
X-ray energies. 

Concerning the main point of this study, namely the 
location of the 7-ray emission site in 3C 279, the equipar- 
tition model makes a rather definite reply for the Epochs 
studied. Provided that one accepts the premise of the 
equipartition model, that the Ce = 1 solution is preferred 
energetically, and assuming t^ar — 10^ s (which can be 
improved), then except for Epoch A, where the emis- 
sion could have taken place inside the BLR, the low Lya 
density implies a location outside the BLR, in the range 
~ 0.1 " 0.5 pc from the central supermassive black hole. 

We also find that for powerful blazars like 3C 279, pro- 
tons satisfy the necessary but not sufficient Hillas condi- 
tion for acceleration above 10^° eV. Though FSRQs are 
not likely to explain the origin of all UHECRs, UHECRs 
in FSRQs can make a VHE 7-ray emission component 
such as has been detected in sources like 3C 279, PKS 
1222-^216, and possibly PKS 1510-089. Such a feature 
would be difficult to explain in the framework of the lep- 
tonic equipartition model presented here, but the uncer- 
tainty related to the variabilty timescale and equiparti- 
tion relation makes it hard to draw definite conclusions. 
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